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Gear box dynamics is characterised by a periodically changing stiffness and a backlash which 
can lead to a loss of the contact between the teeth. Due to backlash, the gear system has piece- 
wise linear stiffness characteristics and, in consequence, can vibrate regularly or chaoticaly 
depending on the system parameters and the initial conditions. We examine the possibility 
of a nonfeedback system control by introducing a weak resonant excitation term and through 
adding an additional degree of freedom to account for shaft flexibility on one side of the gear- 
box. We shall show that by correctly choosing the coupling values the system vibrations may 
be controlled. 

1.1 Introduction 

Gear box dynamics is based on a periodically changing meshing stiffness complemented by 
a nonlinear effect of backlash between the teeth. Recently, their regular and chaotic vibra- 
tions have been predicted theoretically and examined experimentally fT"2"3l|4]|5l|6|[7j. The 
theoretical description of this phenomenon has been based mainly on single degree of free- 
dom models O] |2| E |4j |5j |6| |8l or multi degree models neglecting backlash l9l 1101 . This 
paper examines the possibility of taming chaotic vibrations and reducing the amplitude by 
nonfeedback control methods. Firstly, we will introduce a weak resonant excitation through 
an additional small external excitation term (torque) with a different phase II 1II12I . Secondly, 
we will examine the effect of an additional degree of freedom to account for shaft flexibility 
on one side of the gearbox |7 1. This may also be regarded as a single mode approximation to 
the torsional system dynamics, or alternatively as a simple model for a vibration neutraliser or 
absorber installed in one of the gears. 

1.2 One Stage Gear Model 

We start our analysis with modelling the relative vibrations of a single stage transmission gear 
f6J . Figure fLTI shows a schematic picture of the physical system. The gear wheels are shown 
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Figure 1.1: The physical model of a one stage transmission gear system. 



with moments of inertia Ii and I2 and are coupled by the stiffness and damping of the teeth 
mesh, represented by kz and cz- 

The equations of motion of the system may be written in terms of the two degrees of 
freedom, 4'2 that represent the rotational angles of the gear wheels. These angles are those 
that remain after the steady rotation of the system is removed 0. Thus, if the backlash is 
initially neglected, we have 



-fiV^i + - ^2V'2) + cz(ri?/>i - r2?/'2) ri = Mi, 



(1.1) 



where ri and r2 are the radii of the gear wheels and the overdot represents differentiation with 
respect to time t. 

In fact, it is possible to reduce the above two equations to one using the new relative dis- 
placement coordinate x — riipi — r2ip2- That coordinate represents the relative displacement 
of the gear wheels at the teeth. The equation of motion for x is obtained by subtracting r2//2 
times first equation from ri/Ii times the second one (Eqs. II. 1> . Thus, in non-dimensional 
form, the equation of motion can be rewritten as 

d2 2C d Hr)g{x,7j) B{t) Bp + B,cos{r + 6) 

— 2 — ^ — 2~ ^ 2 ' (^•2-' 

(XT Lu dr cj^ a; 
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Figure 1.2: The meshing stiffness fc(r). 

where the parameters are easily derived from Eqs. (Il.lll.2> . as 

T — Lot, 

2C = cz[rl/h+rl/h], (1.3) 
B{t) = nMi/h+r2M2/l2- 

k{T), g{x,ri) and B{t) are defined as in reference |6|. Note that the backlash and time 
dependent meshing stiffness have been included by rewriting the meshing force as [6J 

kz [rllh + rl/h] X = fc(r).g(x, ry). (1.4) 

Equations M3\ and MA\ have assumed that the moments on the gear system are composed of 
a sinusoidal moment at frequency lo with a constant offset. 

The meshing stiffness is periodic and the backlash is described by a piecewise linear func- 
tion. Figure [O] shows a typical time dependent mesh stiffness variation /c(t) |6| and the 
backlash is modelled for a clearance 77 as 
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1.3 Vibrations of a Gear System in Presence of a Weak 
Resonance Term 

Now we are going to examine the gear system described by Eqs. ( Il.lll.5t subjected to to an 
additional external periodic excitation with a relatively small value D 

d2 2C d k{T)g{x,ri) Bo + Bicos{t + &) + Dcos{t + &') 

T^^^ — ':r^~^ 2 — = 2 ' ^^■°> 

dr^ LO (IT LO LO'^ 
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Figure 1.3: The amplitude of vibration (a) and tiie maximal Lyapunov exponent Ai (b) versus frequency 
uj obtained by simulation. Tlie initial conditions were assumed to be a;o = —2.0 and vq — xo = —0.5 
for small uj (uo = 0.1) and for each new lu (uJi+i) calculations were performed for 400 excitation cycles. 
Each time new initial conditions (a::o, fo) were chosen as the last pair of values of {x, v) for previous 
uji ({xo,vo)\new = {x,v)\oid)- Curves 'V and '2' correspond to cases with D = and D — 2, 
respectively. 



where the phase angle O' ^ 8. Such inclusion has been discussed earlier in the literature as 
a suitable method of chaos control for Josephson junction circuits fl 11 and Froude pendulum 
motion 1 12|. In those papers the authors claimed that the resonant term, if used adequately, 
can tame or include chaotic motion. We have performed numerical simulations of Eqs. il.2\ 
and jl.6> to highlight the effect of the addition of resonant term D on the basic equation of 
motion (Eq. II. 2^ In particular we were interested in changes in vibration amplitude. The 
type of motion was analysed using Lyapunov exponents obtained by the algorithm of Wolf et 
al. 1131 . System parameters have been used as in f6l: w = 1.5, Bq = 1, Si = 4, 6 = 0, 
C = 0.08, ?7 = 10 while the additional excitation term ll.6l was introduced as £> = 2, 9' = 0.1. 

Figure [T3t shows the calculated vibration amplitude of the relative motions of the gear 
wheels. A, defined as 
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Figure 1.4: The Poincare maps calculated for ten randomly chosen initial conditions for u = 1.^ (a) 
and ui — 1.6 the gray colour corresponds to simulations with respect to (Eq. ll.2> while the black one to 
(Ea. ll.6> . respectively. 



The maximal Lyapunov exponent Ai versus frequency cj, obtained by simulation, is plotted in 
Fig. 11.3b . Note, the initial conditions were assumed to be xq = —2.0 and vq — xq — —0.5 
for small starting oj (oj = 0.1) and for each new uJi+i calculations were performed for 400 
excitation cycles with where new initial conditions were the last pair of values of {x, v) for 
the previous frequency uji. Curves marked by '1' correspond to solutions of Eq. 11.21 while 
those marked by '2' correspond to Eq. 11.61 One can easily note that the amplitude of the 
modified system is slightly larger but simultaneously the system behaves more regularly. In 
most of cases where the original system was in a chaotic state (Ai > 1 for the curve '1', 
Fig. \1.3b ) it vibrates regularly in the presence resonant term (Ai < 1 for the curve '2' in 
Fig. 11.3b ). Unfortunately the taming of the chaotic motion was accompanied by additional 
features visible in Fig. II. 3b . Clearly, the curve '2' describing the amplitude of motion shows a 
number of jumps, i.e. for oj « 1.36 and 1.70. These jumps may be associated with transitions 
to other solutions (attractors) with changing oj. This effect is typical for many nonlinear 
systems but it seems to be more transparent for D ^ 0. To explore further this effect we 
have also calculated Poincare maps for many initial conditions chosen randomly. The results 
for frequencies uj — 1.5 and 1.6 are shown in Fig. 11.41 In this figure, one can see that 
for UJ ~ 1.5 attractors (in gray colour) for D — and (in black colour) for D ^ are 
similar Here, numbers '1' and '2' denote the regular and chaotic attractors, respectively. But 
if we move w to a larger value (uj= 1.6) the regular attractor for D = Q disappears. The 
attractor remains for D ^ 0, but interestingly another chaotic attractor emerges. For some 
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other frequencies additional attractors also exist, which supports the proposed explanation of 
the jumping phenomenon. 



I. 4 Vibrations of a Gear System with a Flexible Shaft 

Here we examine the effect of adding an additional degree of freedom to the gear system (Fig. 

II. 5> by using a flexible shaft on one side of the gearbox |7 1. 

Now the equations of motion of the system may be written in terms of the three degrees of 
freedom, t/ii, ?/'2 and T^a, related to the rotational angles of the gear wheels and disk, and have 
the following form 



l2i^2 - kziriipi - r2ip2) + cz{riijji - r2ijj2) ^2 - fcs('03 - V'2) = -Af2, 
h^s + ^5(^-3 - V'2) = -M3. (1.8) 
As in the previous case (Eqs. II. 1111. 2> the above set of equations may be decoupled by 




Figure 1.5: The physical model of the gear model with an additional degree of freedom introduced by 
a disk fixed by a flexible shaft to one of gear wheels. 
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Figure 1.6: The amplitude of vibration (a) and the maximal Lyapunov exponent Ai (b) versus frequency 
ui obtained by simulations. The initial conditions were assumed to be xio = —2.0, «io = iio ~ —0.5 
and X20 = V20 = X20 = for small uj (uj — 0.1) and for each new uji+i calculations were performed 
for 400 excitation cycles where the new initial conditions were the last pair of values of (ai, vi, a;2, W2) 
for the previous u>i. Curves 'r,'2','3' were obtained for fcs = 0, 0.1 and 1.0, respectively. Four nonzero 
Lyapunov exponents A, (i = 1, 2, 3, 4) versus frequency u) are shown for ks = 0.1 (c) and ks = 1.0 
(d). 
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using the new relative displacement coordinates xi — ri-ipi — ^21^2 and X2 — r2{ip3 — ^'2)- 
Thus we reduce the number of equations from three to two 



d^ 2C d k{T)g{xi,ri) l3iksX2 

-r^xi + — 3-2^1 H n 1 5 — 

ar'' oj dr lo^ lo'^ 

Bp + Bicosjujt + 9) 

d^ I32ksx2 , 2/?3C d /?3fc(r)g(a;i,??) 

3-^2:2 + 5 1 3-2;! H ^ = 0, 

dr^ w dr 

where most of parameters have already been defined in Eqs. I1.3l and [l~4l The rest of them are 
easily derived from Eqs. J1.8I1.9> . as 

Pi = I//2, 

p2 - I//2 + I//3, 

P3 = {rl/l2)/[rl/Ii+rl/l2], (1.10) 
= r2M2/l2~r2Ahlh. 



The simulation results of the above system (Eq. 11.91 Fig. II. 5> are presented in Fig. I1.6l and 
11.71 First of all we have repeated the calculations of the amplitude and Lyapunov exponents. 
The results for the amplitude A with various couplings (fcs) are summarized in Fig. II. 6h . 
Curves 'r,'2','3' correspond to for ks = 0, 0.1 and 1.0, respectively. One can see that a 
suitable change to ks may lead to the lowering of the amplitude. This effect is more efficient 
in the case of ks ~ 1.0. Comparing the corresponding values of the maximal Lyapunov 
exponents one can investigate the chaotic nature of the vibrations (Fig. 11.6b ). Namely, the 
case ks = 0.1 (the curve '2') does not reduce the original chaoticity in region around a; « 1.0 
and io G [1.24, 1.62] (the curve '1') while ks = 1.0 is large enough to tame the chaos in these 
regions as well as to make the vibration amplitude smaller. 

It should be noted that for ks ^ Q there are four nonzero exponents to be examined. The 
positive value of the largest one Ai, is presented in Fig. 11.6b . It detects frequency regions of 
chaotic vibrations. But it is also interesting to see what has happened to the other exponents, 
especially for our multidimensional system. They were also calculated and are plotted for 
comparison in Figs. I1.6t -d. For a small coupling stiffness, ks = 0.1, we have the surprising 
result that most of the chaotic regions are characterized by two positive exponents (Fig. II. 6t ). 
which is a signal that the system is truly hyperchaotic II14II15|[T01 . In this case two initially 
close trajectories escape exponentially in two different directions. The results of calculations 
with larger ks {ks = 1.0) are different (Fig. II. 6> . where regions of chaotic motion with 
only one positive Lyapunov exponent were found, which signals typical chaotic behaviour. 
To clarify this point we plot the phase portrait and Poincare maps for the chosen frequency 
uj = 1.5. Figure [l~7k shows the portrait plane and the Poincare map of chaotic attractor for 
ks = while in Fig. 11.7b . for ks — 1.0, the corresponding attractor is regular, synchronized 
with the excitation frequency u. Figure s [l~7t and d show the phase portrait and the Poincare 
map of the hyperchaotic attractor. 
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This paper has examined the effect of adding an additional small resonance excitation term 
or an additional degree of freedom to a simple model of gear vibration (Figs. ll.ll and [L2t . 
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We have shown that the resonance term can successively reduce chaoticity in the dynamical 
system (Fig. 11.3b ). although the vibration amplitude increases to a larger value than that for 
the basic system. We also noted several sudden jumps in the vibration amplitude value (Fig. 
II. 3h ) which could destroy a real gearbox. The jumps are presumably caused by the creation 
of a larger number of attractors in the presence of the resonance term (Fig. II .4> . The extra 
degree of freedom (Fig. I1.5> . which may represent a flexible shaft or a vibration neutraliser, 
also has a considerable effect on the dynamics. We noted that the suitable use of the coupling 
value, shaft stiffness ks, can lead to regular attractors {ks = 1.0) and, simultaneously, to a 
reduction in the vibration amplitude (Fig. I1.6h -b). Interestingly we have found that for a small 
value of the coupling stiffness, ks = 0.1, two of the four Lyapunov exponents were positive, 
leading to hyperchaos phenomenon. The study of the hyperchaotic attactor will be reported in 
a separate article. Thus the coupling ks is a proper bifurcation parameter giving a wide range 
of system behaviour. The correct choice of this stiffness value can also be used to control the 
system vibrations. 
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